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Abstract. We report our detailed calculation of the infrared spectra of GaAs and GaP from first principles, and similar 
ongoing efforts for diamond-type spectra, which have recently been calculated. 



LATTICE DYNAMICS 

Over the last decade or so, density-functional perturba- 
tion theory (DFPT) and application of the 2n+l theorem 
have revolutionized the theoretical study of lattice dy- 
namics [ 1 ]. Simultaneously, sophisticated frozen-phonon 
algorithms have been very successful in calculating com- 
plete phonon dispersion relations|2]. These advances al- 
low detailed knowledge of the vibrational continuum, to 
complement the realistic calculation of electronic band 
structure that has been possible for some time. 

Going beyond the harmonic approximation, with 
DFPT it has been possible to calculate anharmonic 
optical effects, such as widths, shapes, and shifts of 
Raman modes |3], and two-phonon absorption spectra 
in silicon and germanium |4j] . Higher-order anharmonic 
effects can be treated by combining DFPT and frozen- 
phonon methods, and developing both methods may 
prove important to widening the scope of multiphonon 
interactions and their effects which can be realistically 
studied from theory. 

Our focus has been on the frozen-phonon calculation 
of phonon-phonon and photon-phonon interactions. We 
have calculated the absorption spectra for GaAs and GaP 
from the far infrared through twice the frequency of the 
dispersion oscillator, or the zero-momentum transverse- 
optical phonon. The results demonstrate temperature de- 
pendence of the response in the far infrared, the asym- 
metric absorption profile of the dispersion oscillator and 
rich spectral structure caused by infrared light coupling 
to two-phonon states. The correspondence with experi- 
mental spectra are shown in Figure 1 . We have also been 
able to relate prominent features in the two-phonon ab- 
sorption spectra to very simple kinematical parameters 
of the primitive cell, the mass ratio of the basis ions in 
particular. This work combines calculations of Raman 
widths-as has been accomplished with DFPT 0|-with 
the modern theory of polarization 1 5] to calculate absorp- 
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FIGURE 1. The one-phonon (line) and zero-momentum 
two-phonon (+) densities of states are plotted in top panels. In 
bottom panels, the calculated (line) and measured (crosses) ab- 
sorption coefficients are plotted at 300 K for GaAs, and < 10 K 
for GaP. To point out the far-infrared temperature dependence, 
the theoretical spectra, in the far infrared, are plotted at low 
(GaAs) and room (GaP) temperatures (light dots). Plotted the- 
oretical spectral features have been moved to 0.98 of their cal- 
culated position to facilitate comparison with the experimental 
spectra. The data are from Ref. 11(11 . 



tion spectra completely ab initio. 

Currently we are working on a frozen two-phonon cal- 
culation of the weak infrared absorption in diamond-type 
materials. Preliminary calculations of this nature have 
been reported ((J, and at this conference, a thorough cal- 
culation of the spectra for Si and Ge with DFPT is pre- 
sented [4]. A point worth emphasizing is that, with our 
two-phonon calculation, the adiabatic evolution gives no 
change in the ionic moment with respect to the supercell 
origin; that is, the ionic contribution to the "current" is 
zero, because the the ionic displacements are given by 
standing waves. 

All the work presented here is within the frozen- 
phonon approach. In particular, the dynamical matrix and 
anharmonic matrix elements are calculated with wave- 
commensurate supercells, where a phonon's periodicity 



is mapped to a supercell |7]. Born-von Karman force 
matrices can then be calculated as finite differences of 
Hellmann-Feynman forces. The technique has the advan- 
tage that force matrices are sampled in reciprocal space, 
and then interpolated, as is also the case in applications 
of DFPT This allows implementation of the LO-TO 
splitting according to long-established formalisms |8]. It 
also avoids any aliasing. 



ZINCBLENDE CRYSTALS 

Each of the above infrared-active processes is related 
to the interaction between zero-momentum two-phonon 
states and some macroscopic coordinate. In the case of 
GaAs and GaP, within the present theory, it is only the 
dispersion oscillator which is directly coupled to the 
macroscopic field, and the two-phonon states acquire 
oscillator strength through their anharmonic coupling to 
the dispersion oscillator. The relevant matrix element 

for this coupling is M aj3k = £„/,-,-/ ^^ ^e^f. The 
derivative of the dynamical matrix, DT.J (k), is taken with 
respect to the sublattice relative coordinate, and phonon 
polarizations are given on the right. Phonon branches 
are denoted a and j3, Cartesian directions i and i', and 
basis ions t and %' . The matrix element above is actually 
the cubic term in the Born-Oppenheimer potential which 
gives interaction between three normal modes of wave 
vectors k, — k, and 0. 

The two-phonon states, in addition to imprinting two- 
phonon character upon the spectra, also act as decay 
channels, giving the dispersion oscillator absorption fea- 
ture finite linewidth and asymmetric shape. Hence, the 
anharmonic theory of infrared spectra in polar materi- 
als accounts for a broad class of infrared optical prop- 
erties, from sensitive temperature dependence of the far- 
infrared absorption, to the profile of the reststrahlen and 
the two-phonon features. The dielectric function is given 
byH: 
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The numerator on the right side contains the Born ef- 
fective charges, and the denominator contains the ionic 
masses. The frequency of the dispersion oscillator and 
its self-energy appear in the quantity in parentheses. The 
self-energy is temperature-dependent, and is calculated 
via diagrams |9] and the matrix elements defined above. 
The results for GaAs and GaP are plotted in Figure 1, 
along with the experimental spectra. 



DIAMOND-TYPE CRYSTALS 

Since the Born effective charges, which express the inte- 
grated oscillator strength in Eq.^ are zero m a diamond- 
type material, the absorption spectra follow from a dif- 
ferent type of interaction. In fact, the electric field at 
infrared frequencies is a macroscopic coordinate inde- 
pendent of the zero-momentum optical phonon coordi- 
nates, and the key quantity for the spectral calculation 
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is, — — , where E is the electric field. Such calcula- 
tions were recently performed with DFPT and methods 
for incorporating macroscopic fields into electronic cal- 
culations |4]. 

Another approach toward a calculation of the spec- 
tra is to evaluate the polarization change due to a 
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two-phonon displacement: AP, = j t , k 3 J k ^ M ka M -k/3> 
with the modern theory of polarization. Initial work 
in this direction has been reported @. Using wave- 
commensurate supercells, ionic displacements represent 
standing waves, and are appropriate to the momentum- 
conserving two-phonon geometry. For the diamond-type 
structure, there are six branches, and hence the polariza- 
tion is a 6 x 6 matrix at any particular wave vector, rep- 
resenting the possible combinations of phonon branches, 
making the desired quantity a third-rank tensor. We are 
attempting to calculate this tensor in a Cartesian basis. 
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